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We propose a physical implementation of time and spatial parity transformations, as well as 
Galilean boosts, in a trapped-ion quantum simulator. By embedding the simulated model into an 
enlarged simulating Hilbert space, these fundamental symmetry operations can be fully realized and 
measured with ion traps. We illustrate our proposal with analytical and numerical techniques of 
prototypical examples with state-of-the-art trapped-ion platforms. These results pave the way for 
the realization of time and spatial parity transformations in other models and quantum platforms. 
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I. INTRODUCTION 


In the last decade, observing quantum phenomena that 
are difficult or even impossible to detect in the laboratory 
has been possible through the concept of quantum simu¬ 
lation [L . Originally an idea of Richard Feynman 0, it 
is based on implementing a complex quantum dynamics 
on a controllable quantum system. Many proposals and 
experiments on quantum simulations in different control¬ 
lable platforms such as trapped ions 00, superconduct¬ 
ing circuits 00, ultracold gases 0ES, quantum pho¬ 
tonics systems mm, and optical lattices mm, have 
been performed and have led to a deeper understanding 
of a wide variety of phenomena. 

Up to now, proposed models and experimental real¬ 
izations of quantum simulations with trapped ions have 
been realized in spin models OMEU, quantum field the¬ 
ories Q0, quantum phase transitions m, many-body 
systems [20H22I . fermionic and bosonic interactions [2. 'SI . 
relativistic quantum physics, including Dirac equation 
Zitterbewegung , [24] and its realization in the labora¬ 
tory f2Sj, Klein paradox 03 HZ], and interacting Dirac 
particles [28] . among others. Recently, an implemen¬ 
tation of the Majorana equation and unphysical quan¬ 
tum operations was proposed [29] . and experimentally 
realized 01 HP- In addition, U. Alvarez-Rodriguez et 
al. )55] have developed a mathematical formulation of 
an enlarged space or embedding space to perform linear 
transformations between space-time coordinates in a gen¬ 
eral quantum simulator. However, the implementation of 
these concepts in a trapped-ion simulator, including par¬ 
ity V operations, has not yet been analyzed. 

In this Letter, we propose the realization of time and 
spatial parity operations, as well as Galilean boosts, in a 
trapped-ion quantum simulator. We perform analytical 
and numerical calculations in paradigmatic examples to 
illustrate our protocol, which is based on state-of-the-art 
trapped-ion technologies. We show that this proposal, in¬ 
cluding state initialization, dynamics, and measurement, 
can be efficiently implemented in current experiments. 
While Ref. [32] focuses on the underlying mathematical 
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FIG. 1. (color online) Scheme of the proposed experiment for 
the implementation of time, spatial, and Galilean transforma¬ 
tions. Two ions are needed for the time parity and Galilean 
boost, while a single ion suffices for the spatial parity. 


properties behind the theoretical protocol for Galilean 
transformations, here we propose a realistic implemen¬ 
tation in trapped-ion systems. Moreover, we develop 
a toolbox for trapped ions that can be used to imple¬ 
ment reference frame transformations for any given dy¬ 
namical equations. New interesting types of simulations 
can emerge by adding the reference frame transforma¬ 
tion to the toolbox of possible operations. This work 
significantly advances the field of quantum simulations 
of unphysical operations and establishes a path for im¬ 
plementing time and spatial parities in quantum optics 
systems. As a further scope, these results may allow us 
to enhance our capabilities when studying many-body 
interacting systems and their symmetries. 

The formalism introduced in Ref. [32] allows us to 
implement the quantum simulation of reference frame 
transformations in the lab. This symmetry transforma¬ 
tion is described by a linear relation between the initial 
(t, x) and the final {t',x') coordinates, x[ = a ij x j, 
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i,j = 0,1. The spinor in the enlarged space is de¬ 
fined as ^(x,t) = (i/) e , ip°) T , where the even and odd 
part of any wave function can be expressed as ip e ’° = 
± rp(x',t')\. Therefore, the dynamical informa¬ 
tion of ip(x, t) and ip(x', t') is encoded in the evolution of 
d'. Moreover, through a judicious choice of measurement 
observables, one can perform a reference frame transfor¬ 
mation via a local cr z operator, or even observe spacetime 
correlation functions between different reference frames. 
Throughout this paper we consider the evolution in the 
simulated Hilbert space as given by the Schrodinger equa¬ 
tion idti/) = —icd x ip , where c is a simulated speed of light 
and we fix h = 1. The corresponding equation for tl/ in 
the embedding space, for arbitrary Galilean transforma¬ 
tions, may be written as 

idt^ = -i(5il + a 2 o x )d x ^, (1) 

where di t 2 = [c(an ± aoo) T Q'io]/(2q'h). We explain 
now how to use this representation for the implementa¬ 
tion in a trapped-ion system. In the Lamb-Dicke regime, 
r\\J ((° + cd) 2 ) •C 1, the Hamiltonian describing the in¬ 
teraction between an ion and a laser driving is (3. 

U{t) = fl of 7+[l + in(ae~ ivt + e ivt )]e i <*- st '> + H.c., 

where 5 = ui — ojq is the laser detuning, r/ = k^l/2mv is 
the Lamb-Dicke parameter [Jj, & is the wave number of 
the external field, m is the mass of ion, v is the frequency 
of a static potential harmonic oscillator, fl 0 is the cou¬ 
pling strength, a and od are the annihilation and creation 
operators of any suitable vibrational mode of the ion 
string, that we choose to be the center of mass motional 
mode, and 4> is the field phase. In spin-1/2 language, 
cr+ = \ e )(g\ = {a x +io v )/2, a~ = \g){e\ = (o x - io y )/2. 

When S = 0, a carrier resonance can be obtained with 
Hamiltonian H c = Q(a + e 1 ^ + cr~e -i ^). A red-sideband, 
also known as Jaynes-Cummings (JC) interaction, is re¬ 
alised in the case of S = —v. This Hamiltonian is written 
as H r = Q,g{ao + e 1 ^ 1 ' + a^o~e~ l ^ r ). Respectively, when 
S = u, a blue-sideband, anti-Jaynes-Cummings (AJC) 
interaction can be achieved, and its Hamiltonian can be 
expressed as Hj, = Clr](a^o + e l,l>b + aa~e~ l(l>b ). We illus¬ 
trate now how to apply these techniques to generate the 
time and spatial parity transformations in trapped ions. 


II. TIME PARITY TRANSFORMATION 

As a first example, we show how to use two trapped 
ions to simulate a time parity transformation, (t, x) —> 
(~t,x), (a 0 o,a 0 i,aio,Q:ii) = (-1,0,0,1). Here, we 
choose the Hamiltonian in the simulated space as a time- 
independent one, H = H e = cp and the momentum 
p > 0, which describes a massless Dirac Hamiltonian 
without the internal degree of freedom. The correspond¬ 
ing one-dimensional Schrodinger equation in the enlarged 
space can be expressed as 

( 2 ) 


The Hamiltonian in the enlarged space, H = erf ® H = 
ofcp, where erf is the Pauli operator acting on ion 1, 
which can be realized by implementing a blue- and red- 
sideband simultaneously [3J [M], 

T-L = r/Q(<Ti a^ 1 e^h+crf ae- l ‘t >b )+ri&((j+ae l ' l,r +<7f 

(3) 

with proper phases for blue- and red-sideband (f>b = 7t/2, 
= —7t/2. Here, ryD = ^ and i(af — a)/2 = pA with 

A = a/1/2 mv. We depict in Fig. [l] a scheme of the 
experimental setup with two ions interacting with lasers. 
The initial state in the enlarged space is given as, 

f = 0) = ® ^{x, t = 0), (4) 


where ip(x,t = 0) can be described as a Gaussian 
wave packet, ip(x,t = 0) = ipo{x,t = 0)e ip ° x = 
(^/\/^nK)~ 1 e~ x2 / 4A2 e lp ° x . In a trapped-ion setup, this 
can be achieved by cooling the motional mode to the 
ground state, which is a Gaussian, and displacing it by 
simultaneous red and blue sidebands with the Hamilto¬ 
nian poxaf , where x = A (a + cd), with different phases 
in the a and cd operators as compared to the simulating 
Hamiltonian T-L. This will allow one to achieve an average 
Po momentum, using the auxiliary second ion initialized 
in an eigenstate of the of operator. After applying the 
evolution propagator exp(— iHt), we can evolve the state 
for any time. The solution reads 
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Furthermore, the quantum states in the simulated spaces 
are obtained reversing the initial mapping, 

ip{x,t) = (1, l)4 / (a:,t) = e~ ( 4A 2 ' e -wo(ct-x) ^ 

V \[2ttA 

i p(x,—t) = (1, l)o z ty(x, t) = e e ip 0 (ct+x)' 

V y/2wA 


( 6 ) 


We plot in Fig. [2](a) the initial wavepacket in the simu¬ 
lated space, and in Fig. |2](b) the evolved and time-parity- 
transformed wavepackets in Eq. We calculate now 
the position average values in the simulated space for the 
different reference frames and their correlation, 


{&}il>(x,t) {fc)il>(x,—t) ^d t) 0* 

( 7 ) 

The spacetime correlations may be measured in dif¬ 
ferent ways with current ion technology. The one we 
introduce here extends the physical principle employed 
in (25) . The measurement is performed upon the of 
observable of the first ion, associated with the enlarged 


idt'S’ = ofcp' F. 
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space degree of freedom, via fluorescence detection. To 
achieve this, a state-dependent displacement operator 
U = exp(— ikxaf/2) is applied to the internal state 
of this ion and the joint mode, with A = U^afU = 
cos (kx)af + sin (kx)a\ [J5 . In order to detect, e.g., the 
spacetime correlation (x)^^ x ^)^{ x -t) in the simulated 
space, one should measure (x(af + t)a i)^( x ,t) follow¬ 
ing Eq. |6]). Here, the operator (erf + l)crf acts on the 
qubit degree of freedom of the enlarged space, and the 
measurement can be decomposed into two parts, one for 
each summand, (xa and — i(xa\)^i xt y These 
two measurements can be obtained from the derivative of 
the A observable with respect to k in the limit k(x) <C 1, 
in which dk(A) ss (xaf), with a local rotation in the first 
case to change erf into erf. Moreover, computing (A) for 
sizable fc, for the cases of initial erf and cr\ eigenstates 
in the internal state, allows one to obtain cos (kx) and 
sin (kx). Via Fourier transform, we can obtain the posi¬ 
tion wavepacket probability distribution. The previous 
procedure enables us, among other things, to compute 
spacetime correlation functions without full tomography, 
which can reduce significantly the required resources. 
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FIG. 2. (color online) (a) Probability distribution |i/;(a:)| 2 
at time t = 0. Probability distribution \%jj(x)\ 2 at time t = 
2 x 10 3 A with c = 1CF 3 and v/c = 0.8 for time parity (b), 
spatial parity (c), and Galilean Boost (d). Dashed blue lines 
are time evolutions and solid red ones are the corresponding 
Galilean group transformations. 


III. SPATIAL PARITY TRANSFORMATION 


The second case we consider is the simulation of 
a spatial parity transformation, {t,x) -A (t,—x), 
(aoo,aoi,ctio,an) = (1,0,0,—1). The initial state 
in the simulated space coincides with the one of time 
parity, ip(x,t = 0) = (\/V^nA)~ 1 e~ x2 / 4A2 e lp ° x . So 
does the Hamiltonian in the simulated space, H = 
H e = cp. It is obvious that ij)(—x,t = 0) = 
(\J\ / '2nA)~ 1 e~ x / 4A e~ lpoX . The Hamiltonian in the en¬ 
larged space H goes to erf G> H e , and the initial spinor in 
the enlarged space can be expressed as 


V(x,Q) 


1 _ 

- C 4A= 

2 v v^brA 


e ip 0 x _|_ e ~ip 0 x \ 
e ip 0 x _ e ~ip 0 x J 


( 8 ) 


This can be achieved by initializing the internal state 
associated with the enlarged space in the (1, 0) T state, 
cooling the motional mode to the ground state, which is 
a Gaussian, and performing a conditional displacement 
of the motional state with the Hamiltonian p 0 xaf. We 
point out that for spatial parity one ion suffices for both 
state initialization and simulation. 

Accordingly, the state in the enlarged space consider¬ 
ing time-evolution results, 


^(M) = 


2\A/^A 


(9) 


(< 4A^ e ~iP0(ct~x) _|_ g- <C 4^2 ) e ~ip 0 (ct+X) 
(Ct 4A V p-ipn(ct-x) „~ (ot AV „—ipn(ct+x) 


As in the previous case, in order to recover the wavefunc- 
tion in each of the simulated spaces we add or subtract 


the even and odd components of the spinor 11/. 

ip(x,t) = . 1 ( 10 ) 

V y/2rrA 

*/>(-x,t) = 1 ! e -i^r g e- ip °( ct + a! ). (ii) 

V y/2n A 


We plot in Fig. |2jc) the evolved and spatial-parity- 
transformed wavepackets in Eqs. (10) and We also 
obtain the expectation values in each of the correlations. 


(x}ip(x,t) ct 7 x,t) 


= (T(a;,t)|(Tfa;| v I , (x,t)) — i('I>(x,t)\a\x\'ff(x,t)} 

= (\H(:r,t)|crfa;|\E ’(x,t)) — i($(x, t)\(T*x\$(x, t)) 

= —2ipoA 2 e~ j* 2 e _2p ° A , (12) 

where <f>(a;,f) = e -171 ’ 0 ’ 1 ^^(x, t). These can be measured 
as in the previous example. We point out that the spatial 
parity case can be distinguished from the time parity case 
through the spacetime correlation, which is different in 
both cases. The computation of these correlations with 
our method makes full state tomography unnecessary. 


IV. GALILEAN BOOST 

Here, we propose the simulation of a reference frame 
change associated with a Galilean boost, (t, x) -A (t, x — 
vt), (aoo, aoi,aio, a n) = (1, 0, —v, 1), where v is the rel¬ 
ative velocity [33]. Here, we consider the previous Hamil¬ 
tonian and initial state in the simulated space. The cor¬ 
responding Hamiltonian in the simulating enlarged space 
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reads 


H(t) = (c+^jtp- l -a x p. (13) 


Hi (t) 


H.2 (t) 


Moreover, we can calculate the initial spinor state and 
compute the time evolution. The expression for the quan¬ 
tum states in the simulated spaces reads, 


= —-- e ~ lv ^ ct ~ x \ (14) 

V \/ 27 r A 

ip{x - vt,t) = , 1 e~ (ot ~^ Vt)2 e- lMct ~ x+vt \ (15) 

V \/27rA 


We plot in Fig. [2](d) the evolution of the wavepackets 
with and without Galilean boost in Eqs. (14) and (15). 
Moreover, the expectation values for position x are 


{&}p(x,t) — cf, vt,t) — (c -(- v)t 1 (Id) 

{x)i,(x,t),i>(x-vt,t) = ^t{2c + v)e~T& e~ ipatv . (17) 


For the trapped-ion simulation the initialization of the 
spinor can be done similarly to the time parity case. For 
the subsequent dynamics we divide the Hamiltonian in 
Eq. (13) into two parts to implement its evolution in the 
trapped-ion system. To realize "Hi in the laboratory, we 
propose to use a second auxiliary ion initialized in an 
eigenstate of a%, 


m\n+) = (c+1) ipi*>i+> = ( c + 1) ofp|*>i+>. 

( 18 ) 

Then, the Hamiltonian can be implemented as 

'Hi = ofp (c+ = iri&iaZi^ - a), (19) 


with rj fii = (c + |)/2A. Moreover, the second term in 
Eq. <fT3l) can be realized as 


%2 = ~\ a iP = (« f - a), (20) 


and with pCl 2 = — , through simultaneous red and blue 

sideband excitations upon the first ion. 


1 
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FIG. 3. (color online) Fidelity F = Tr[p\ipi) (ipi\] of (a) time 
parity, (b) spatial parity, and (c) Galilean boost with SI = 
O.Olu (red, upper), SI = 0.025^ (black, middle) and SI = 0.04^ 
(purple, lower), p denotes the final state after initialization 
and dynamics, evolved with Eq. (21), and \ipi) denotes the 
ideal evolved state in absence of imperfections. We point ou t 
that Cl in the Galilean boost case (c) equals S7i in Eq. (191. 
We consider the parameters of trapped Ca + ion in some of the 
Innsbruck experiments ESI p = 0.06 and realistic decoherence 
rates [H IA = F c = T_ = 3.7 x lCU 7 ^ and = 6.2 x 10 _7 ja 


V. DISCUSSION 


To analyze the robustness of the simulating system, we 
computed the dynamics with a master equation including 
different decoherence sources. We considered unintended 
carrier transitions due to off-resonant coupling, heating 
Tft,, phonon loss r c , dephasing and spontaneous emis¬ 
sion r_, 

p =-i[T-L T , p] + Y h L(a^)p + Y c L(a)p 

+Y^L(a~~)p + Y-L(a-)p, (21) 


where the Lindblad superoperators are L{X)p = 
(2 XpX^ — X { Xp — pX^X)/2. Here, Hr is the trapped- 
ion Hamiltonian corresponding to each of the three 
cases analyzed, namely, time parity, spatial parity, and 
Galilean boost. We include carrier and counterrotating 
sideband terms in the dynamics, i.e., without perform¬ 
ing vibrational rotating-wave approximation. Therefore, 
this master equation accounts for all the significant deco¬ 
herence and error sources present in current trapped-ion 
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FIG. 4. (color online) Fidelity of (a) time parity, (b) spatial 
parity, and (c) Galilean boost with Cl = 0.0W (red, upper), 
Cl = 0.025^ (black, middle) and Cl = 0.04^ (purple, lower),by 
considering trapped Be + ion as an example. We use the pa¬ 
rameters in some of the NIST experiments [4[: Lamb-Dicke 
factor p = 0.3, decoherence rates Th = r c = T_ = 3.7x lO -7 ^ 
and r^, = 6.2 x lO^ 7 ^. 


experiments.We plot in Figs. [3|a)-(c) and Figs. @a)-(c) 
the fidelities of trapped calcium and beryllium ions af¬ 
ter state initialization and evolution with the dynamics 
in Eq. (21) for the cases of time parity, spatial parity, 
and Galilean boost. For the state initialization part, we 
compute the dynamics with an equivalent master equa¬ 
tion for the corrresponding initialization Hamiltonian as 
described for each case (a)-(c) in the text. We consid¬ 
ered for the initialization time poA = r/Clt = 1 in all 
cases. We point out that because of the relatively small 
atomic mass of beryllium, large trap frequency, sizable 
Lamb-Dicke factor and the corresponding large decoher¬ 
ence rates are introduced in the realistic experiments of 
NIST [lj. Our work shows a significant feasibility in dif¬ 
ferent trapped-ion setups. 


VI. CONCLUSIONS 


To summarize, we have proposed the physical imple¬ 
mentation of fundamental symmetry transformations in¬ 
cluding time and spatial parity with trapped ions. The 
formalism permits as well to perform Galilean boosts 
with the same technology. By embedding the simulated 
physical system into an enlarged Hilbert space living in 
the trapped-ion system, the proposed formalism can be 
carried out with current ion-trap setups. Furthermore, 
our work establishes a path for the realization of parity 
transformations in other quantum platforms and many- 
body interacting models. 
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